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Abstract 

The magnetic dipole and electric quadrupole hyperfine constants of Aluminium ( 27 Al) atom are 
computed using the relativistic coupled cluster (CC) and unitary coupled cluster (UCC) methods. 
Effects of electron correlations are investigated using different levels of CC approximations and 
truncation schemes. The ionization potentials, excitation energies, transition probabilities, oscillator 
strengths and nuclear quadrupole moment are computed to assess the accuracy of these schemes. The 
nuclear quadrupole moment obtained from the present CC and UCC calculations in the singles and 
doubles approximations are 142.5 mbarn and 141.5 mbarn respectively. The discrepancies between 
our calculated IPs and EEs and their measured values are better than 0.3%. The other one-electron 
properties reported here are also in excellent agreement with the measurements. 

PACS number(s) : 31.15.Ar, 31.15.Dv, 31.25.Jf, 32.10.Fn 

1 Introduction 

Theoretical studies of properties like hyperfine coupling constants and transition probabilities are strin- 
gent tests of the accuracies of atomic wave functions. The former is sensitive to the nuclear region while 
the latter crucially depends on the wavefunctions at large distances. High precision calculations of these 
properties require a rigorous incorporation of correlation effects and in some cases even relativis- 
tic effects. In particular, the hyperfine coupling constants and transition electric dipole (El) moments 
calculations are relevant to the studies of parity non-conservation (PNC) in atoms as PNC transition 
amplitudes involve both short range electro- weak interaction and El transition moments [2]. 

The relativistic and electron correlation effects can be incorporated in many-electron systems through 
a variety of many-body methods. Among these approaches, the relativistic coupled cluster (RCC) method 
has emerged as one of the most powerful and effective tool for a high precision description of electron 
correlations in many-electron systems p||. Coupled-cluster (CC) is an all-order non-perturbative theory, 
and therefore, the higher order electron correlation effects can be incorporated more efficiently than 
using the order- by-order diagrammatic many-body perturbation theory (MBPT) 0]. The CC method is 
size-extensive [^J , a property which has been found to be crucial for an accurate determination of state 
energies, bond cleavage energies for molecules and related spectroscopic constants. Since the order- by- 
order MBPT expansion terms are directly related to the terms in the CC wavefunction (as the latter is 
an all-order version of the former scheme), the CC results can be improved by adding certain important 
omitted diagrams by computing the corresponding low order MBPT diagrams to all order. 

In this paper, we report our calculations of the magnetic dipole and electric quadrupole hyperfine 
constants (A and B respectively) for the lowest two 2 P^/2 states (3 2 P 3 / 2 and 4 2 P 3 / 2 ) of 27 Al obtained 
using the RCC method. We also present ionization potentials (IPs), transition energies (EEs), transition 
probabilities and oscillator strengths of 27 Al. Effects of electron correlations on these quantities are 
investigated using different levels of CC approximation. We compare atomic properties of 27 Al obtained 
from CC and UCC methods to assess the relative performance and accuracy of these two schemes. 
The UCC and its variants @ |H1 E] were developed almost two decades ago to incorporate higher 
order electron correlation effect systematically. Recently we have applied the relativistic UCC to atomic 
systems for the first time to calculate properties like lifetime of excited states ^H]. To our knowledge, 
no prior UCC calculations are available for 27 Al. 
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The nuclear quadrupole moment (Q) of 27 Al is of interest in several research areas [U]. The elec- 
tric quadrupole hyperfine constant (B) of 27 Al was measured in ionic crystals |12l Il3j and in metallic 
alloys |14; and the value of Q is extracted by combining the calculated electric field gradient (q) at the 
nucleus with the measured value of B. Q is also obtained from studying AlF and AlCl molecules [15] . 
Pernpointner and Visscher ^H] have obtained the value of Q for Al, by studying AlF, AlCl and AlBr 
molecules using fully relativistic CCSD(T) theory. The value of Q is also obtained from the muonic x-ray 
[171 118j and nuclear scattering experiments \W\ . 

In 1976, Rogers et al [201 employed the second order MBPT method to determine the nuclear 
quadrupole moment Q of 27 Al. Later, Sundholm and Olsen j2l] calculated Q for the 2 Pz/2 state of 
Al using the multi-configuration Hartree-Fock (MCHF) approach [22]. Nuclear structure calculations of 
Q have also been carried out |23l 1241 l2"51 12(i| . The discrepancies between the calculated and measured 
values of Q suggest that inclusion of higher order electron correlation effects is necessary to improve 
the existing calculations. Our present work is motivated by this consideration. In this work, we have 
compared our calculated Q value of the with all the available calculated and measured values. 

Section [21 briefly reviews the CC method. Computational details and results are discussed in the 
subsequent sections El and 01 respectively. Finally in the last section we highlight the findings of our work. 



2 Methodology 

Since the coupled cluster methods used in this work are discussed elsewhere El El E3 EE] m details, 
we only outline the essential features of the method here. 

In this work, we employ the straight forward extension of non-relativistic coupled cluster theory to 
the relativistic regime by adopting the no-virtual-pair approximation (NVPA) along with appropriate 
modification of orbital form and potential terms [221 ■ We begin with Dirac-Coulomb Hamiltonian (H) 
which is expressed as 

N N 2 

H = Y, [coh ■ Pi + Pmc 2 + V N { ri )] + — ■ (1) 

■ i - ^ ■ ^7 

i—l i<.j J 

The normal order form of the above Hamiltonian is given by 

H = H N - (0|F|0> = X>'|/|j) {afo} + Ml**) {»M a '°*} . ( 2 ) 

ij i,j,k,l 

where 

(ij\\kl) = (ij\ — \kl)-{ij\ — \lk). (3) 

The valence universal Fock space open-shell coupled cluster method is employed which begins with the 
decomposition of the full many-electron Hilbert space of dimension N into into a reference space Mo of 
dimension M <C N, defined by the projector P, and its orthogonal complement A4q associated with the 
projector Q = 1 — P, A valence universal wave operator ft is then introduced which satisfies 

\9i)=n\*W), i=l,...,M (4) 

where |^^) and are the unperturbed and the exact wave functions of the ith eigenstate of the 
Hamiltonian, respectively. The wave operator ft, which formally represents the mapping of the reference 
space A4q onto the target space M. spanned by the M eigenstates l&i), has the properties 

QP = Q, PQ = P, Q 2 = Q. (5) 

With the aid of the wave operator fl, the Schrodinger equation for the M eigenstates of the Hamiltonian 
correlating with the M-dimensional reference space, i.e., 

fr|* < )=^|* i >, i = l,...,M, (6) 

is transformed into a generalized Bloch equation, 

HflP = flHQP = flPH cS P, (7) 
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where H c g = PHVtP is the effective Hamiltonian. Once Eq. is solved for the wave operator Q, 
the energies Ei, i = 1, . . . , M, are computed by diagonalizing the effective Hamiltonian H e s in the M- 
dimensional reference space Mo- Following Lindgren's formulation of open-shell CC j2j], we express the 
valence universal wave operator ft as 

Q = {exp(a)}, (8) 

and a being the excitation operator and curly brackets denote the normal ordering. 

The operator a has two parts, one corresponds to the core sector and the other to the valence sector. 
In the coupled-cluster singles and double (CCSD) excitation approximation the excitation operator for 
the core sector is given by 

T = T X + T 2 = y, KM «s + \ E K<4<W €1 > (9) 

ap abpq 

and t^l being the amplitude corresponding to single and double excitations respectively. In UCC 
theory the core excitation operator has a unitary form and is represented as T — Th For a single valence 
system the excitation operator the valance sector turns out to be exp(S') = {1 + 5}and 

Sk = Si k + S 2 k = E { a l a k} s l + E { a i a i a bak} sH , (10) 

k^p bpq 

where and denotes the single and double excitation amplitudes for the valance sectors respectively. 
In Eqs. © and il 1 Ofl we denote the core (virtual ) orbitals by a,b,c...(p,q,r...) respectively and v 
corresponds to the valance orbital. In the unitary counterpart of CCSD, i.e. in UCCSD, since the core 
excitation operator also contains a de-excitation part (denoted by T^) it can be shown that for a given 
approximation the UCC theory contains certain higher excitations effects which is not present in the CC 
theory (TH]. 



2.1 Computation of one-electron properties 



We now present the method for computing the matrix-element of sum of one-body operator O = YliLi °* 
that utilizes the structure O = {exp(cr)}. In this approach, the CC-equations are first solved to determine 
the a cluster amplitudes and then the matrix-element of a one-body operator is computed through the 
following relation: 

Q (ttj|Q|tti) (u) 

where \^k) denotes the exact k-th state wave-functions. It can be shown that the substitution of 
the expression for the exact wave-functions l^i) and \^ f) in Ea. ljlljl explicitly cancels out spurious 
disconnected terms from the above expression which reduces to 

(*/|Q|^)c (12] 

where subscript c refers to the 'connected' terms. 



2.2 Magnetic dipole and electric quadrupole hyperfine constants 

The interaction between the various moments of the nucleus and the electrons of an atom are collectively 
referred to as hyperfine interactions 21 . Here we will briefly present and outline of the the magnetic 
dipole (^4), electric quadrupole (B) hyperfine constants and the nuclear quadrupole moment (Q). 
For a state \IJFMp)the magnetic dipole hyperfine constant A is defined as 

A- m m , '• ,||T '" I|J> , as) 

where [ij is the nuclear dipole moment defined in units of Bohr magneton I and J are the total 
angular angular momentum for the nucleus and the electron state respectively and F = I + J with the 
projection Mp. The electric quadrupole hyperfine constant B for the same state is defined as 
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B = 2eQ 



2J(2J-1) 



(2 J +1) (2 J + 2) (2 J + 3) 



1/2 



<J||T( 2 >||J), 



(14) 



where Q denotes the nuclear quadrupole moment. 

The single particle forms (t^) of the operator T^ k \k = 1, 2) are taken from Cheng's paper |3f)j and are 
represented as 



^ (1) =E4 X) = E 



(15) 



and 



^ 2) = E4 2) = E- e ^ 2) (^)- 



(16) 



(k) 



(2fc+l) rfc <? 



• In 



Here a? is the Dirac matrix and Y^ q is the vector spherical harmonics and C, 

Eo. ljTs)) the index j refers to the j-th electron of the atom and e is the magnitude of the electronic 
charge. 



2.3 Electric dipole transition probabilities and oscillator strengths 

The transition probability Aj^i (in sec^ 1 ) and oscillator strength faf (in a.u.) for the electric dipole 
allowed transitions are given by pH] 



.4 



/«-» 



2.0261 x 10 18 



s 



/<-< 



and 



fit = 1.499 x 10- 16 ^A 2 ^ 

9i 



(17) 



(18) 



respectively. Here, A is the wave length in A and gf(gi) = (2 J + 1) is degeneracy of the upper (lower) 
level. The quantity Sf^i is the El line strengths (in atomic units), respectively. The line strengths Sf<-i 
is defined as 



S/,-* = D if x Dfi 
where the electric dipole Dfi matrix elements is given by 



(19) 



D fi = C(f,i) / dr[P f (r)P % {r) + Qj^Q^r . 



with 



C(/,i) = (-1)^/2 ^ J> J ^ v /(2j / + l)(2j i + l). 



(20) 



(21) 



3 Computational Details 

The Fock-space relativistic coupled cluster method is applied to compute the ground and excited state 
energies of Al. The Dirac- Fock equations are first solved for the alkali metal ion M + , which defines the (0- 
hole,0-particle) sector of the Fock space. The ion is then correlated using the closed shell CCSD/LCCSD, 
after which one-electron is added following the Fock-space scheme 

M+(0,0) — > M+(0,1). 

Here LCCSD corresponds to linearized coupled-cluster in singles and doubles. Both the DF and relativis- 
tic CC programs utilize the angular momentum decomposition of the wave-functions and CC equations. 
Using the Jucys- Levinson-Vanagas (JLV) theorem (22], the Goldstone diagrams are expressed as a 
products of angular momentum diagrams and reduced matrix element. This procedure simplifies the 
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computational complexity of the DF and relativistic CC equations. We use the kinetic balance condition 
to avoid the "variational collapse" |33j . 

In the actual computation, the DF ground state and excited state properties of Al are computed using 
the finite basis set expansion method (FBSE) p2| with a large basis set of (32s28p24dl5/) Gaussian 
functions of the form 

F tM {r)=r k -e~ a ^ (22) 
with k = 0, 1, . . . for s,p, . . . type functions, respectively. For the exponents, the even tempering condition 

OLi = aoP 1 - 1 (23) 

is applied. Here, N is the number of basis functions for a specific symmetry. The self-consistent DF 
orbitals are stored on a grid. It is assumed that virtual orbitals with high energies do not contribute 
significantly to properties like IPs. In the CCSD calculations, we therefore truncate the virtual s, p, d and 
/ orbitals above 1000 a.u., 1000 a.u., 500 a.u. and 100 a.u., respectively. Single and double excitations 
from all the core orbitals to valence or virtual orbitals are considered. 



Table 1: Transition energies (in cm 1 ) of Al atom. IP is the ionization potential, EE denotes the 
excitation energies with respect to the 2 P\/2 ground state. 





Dominant 
Configuration 


State 


LCCSD 


CCSD 


UCCSD 


Expt35j 


IP 


[Mg]3 Pl/2 


*Pl/2 


48194.92 


48155.42 


48211.83 


48279.16 


EE 


[Mg]3p 3 /2 

[Mg]4s V2 


2 P 3 /2 

2 s 1/2 


133.94 
24802.66 


114.75 
24937.55 


114.55 
24988.00 


112.04 
25347.69 




[Mg]4p 1/2 
[Mg]4p 3/2 


2 Pl/2 
2 P 3/ 2 


32464.50 
32481.26 


32521.12 
32537.68 


32572.34 
32588.94 


32949.84 
32965.67 



4 Results and Discussions 

Table compares the IP and EE of Al computed using different CC methods with the experiment [35] . 
It can be seen from this table that UCCSD calculations of the IP and EEs are more accurate than the 
CCSD and LCCSD results. Although not well understood, the present as well as some earlier studies 
PQ indicate that the IPs computed using the LCCSD scheme are sometimes in better agreement with 
the experiments than the corresponding CCSD calculations. For instance, the 2 -Pi/2(3pi/2) IP estimated 
using LCCSD method is differs by 84 cm" 1 from the measured value, while the corresponding CC value 
is off by 124 cm -1 . However, it is clear from TableQthat CCSD estimates the EE more accurately than 
the LCCSD scheme for the low lying states. 

In table we present the results of our nuclear quadrupole moment (Q) calculation using different 
CC methods with other calculations |2Ql 12 If and different measurements |151 1171 fl~8l I19j . Pyykko has 
reviewed, calculated and measured Q values for a number of systems [ll]. Comparison of our results with 
the existing data will give an indication of the potential of the CC and UCC methods to provide accurate 
estimate of nuclear properties. It is evident from table[3]that Q calculated using the second order MBPT 
calculations by Rogers et al. J2D] is far outside the experimental limits whereas the value obtained by the 
restricted active space (RAS) based multi-configuration Hartree-Fock (MCHF) [22 calculation is closer 
to the experimental limits. The uncertainty in molecular experiment is less compared to the muonic 
experiments. Although our CC and UCC results are slightly outside the experimental limits, they could 
be of some importance in determining the accurate value of Q from a wide range of values. 

The MBPT (2) and MCHF results clearly indicate that the contributions from non linear terms present 
in CC and UCC theories are non-negligible and this is further supported by the results of our different 
CC calculations of Q. The extremely accurate estimate of Q offered by LCCSD scheme is perhaps a 
bit fortuitous. Nevertheless, the performance of CC, especially the UCC, outshines the MBPT(2) and 
MCHF treatments. Note that the effects of partial triple and quadrupole excitations are present in 
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Table 2: Comparison of the nuclear quadrupole moments Q (in mbarn) of 27 A1 estimated using various 
CC approach with the experiment and with other theoretical calculations. 



Method Q 


t r^r^QTi 


140. 1 




14Z.0 


TTrrcn 


141.0 


IVlJjjr 1 [ZJ ZU 


1 AKfO^ 

Luoyz ) 


Mthf [zzj 


14(J.o(l.UJ 


ivioiecLuai iijA.pi 10 


i^O.O^l -U J 


Molecular Theory |16| 


146.0(4) 


Muonic Exp-UHHHj 


150(6) 


Nuclear Scattering 19 


155(3) 


Nuclear TheoryE3EI] 


134 


Nuclear Theorygl EU 


150.8 


Nuclear TheorygHJ 


138.9 



our UCC calculations. The CC and UCC theories unlike the MCHF method are size-consistent and 
incorporate certain higher order excitations that the MCHF does not at the same level of approximation. 
For example, at the level of single and double (SD) excitations, the CC theory includes not only the 
effect of T 2 but also Xf ; whereas the effect of Xf can be obtained in MCHF only if one considers the 
quadrupole excitations. Also, that calculation is non-relativistic with a relativistic correction added while 
our calculation is fully relativistic. 



Table 3: Magnetic dipole hyperfine (A) matrix elements ( in MHz) of Al atom. 



Method 


A(3p 1/2 ) 


A(3p 3 /2) 


A(4s 1/2 ) 


A(4p 1/2 ) 


A(4p 3/2 ) 


LCCSD 


493.30 


108.39 


414.14 


55.97 


26.26 


CCSD 


498.06 


101.49 


405.94 


58.32 


23.09 


UCCSD 


498.33 


100.98 


407.18 


58.28 


23.12 


Expt [36 


502.0346(5) 


94.27723(10) 









The values of A for the ground and excited states of Al computed using the LCCSD, CCSD and 
UCCSD methods are displayed in Tabled Our calculated values of A agrees well with the experimental 
values [36 for the 2 -Pi/ 2 (3pi/ 2 ) and 2 X3/ 2 (3p3/ 2 ) state. We also present the values of A for some other low 
lying states which could be useful for experimentalists. We have also computed the electric quadrupole 
hyperfine constant (B) for two low lying 2 -P3/ 2 (3j?3/ 2 ) and (4p 3 / 2 ) states using CCSD(UCCSD) theory 
which are respectively 19.49 MHz(19.59 MHz) and 2.85 MHz(2.86 MHz) whereas the experimental value 
of B for the 2 J P 3 / 2 (3p 3/2 ) state is 18.915 MHz [S3EHj- 

In Table 01 we compare the "ipi/2 — > 4s and 3p 3 / 2 3s wave lengths (A), oscillator strengths (fik), 
line strengths (Sik) and transition probabilities (A^) obtained from LCCSD, CCSD and UCC schemes 
with the experiment. Table 0] shows that our computed quantities (A, Sik, and Aik) are in excellent 
agreement with experiment especially those predicted by the UCC scheme. That this scheme provides 
more accurate estimates of IP, EE etc. and is also evident from Figure where the absolute errors (in 
%) in the computed properties are plotted against different CC schemes. 

5 Conclusion 

The relativistic open-shell coupled cluster scheme for direct energy difference calculations and several 
one electron properties is presented and applied to Al. In this work, we investigate the effects of electron 
correlations on the ground and excited state properties using different levels of CC approximations. We 
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Table 4: Wave lengths A (in A), line strengths Sik = k| 2 /ao ( m a - u -)„ transition probabilities A%k (in 
10 8 s _1 ), and oscillator strengths fa (in a.u.) for [Mg]3p —* [Mg]4s transitions of Al atom. 



Method 




[Mg]3p 1/2 


^4s 






[Mg]3p 3/2 


-> 4s 






A 


Sik 


A lk 


fik 


A 


Sik 


Aik 


fik 


LCCSD 


4031.82 


3.379 


0.522 


0.127 


4053.72 


6.763 


1.028 


0.127 


CCSD 


4010.05 


3.292 


0.517 


0.125 


4028.59 


6.634 


1.028 


0.125 


UCCSD 


4001.92 


3.275 


0.517 


0.114 


4020.35 


6.600 


1.029 


0.125 


Expt[39 


3944.01 


2.99 


0.493 


0.115 


3961.52 


6.0 


0.98 


0.115 




Figure 1: Absolute error (in %) of the computed IP, 2 Pi/ 2 (3p 1/2 ) ^ 2 S 1/2 (As) wave lengths (A), magnetic 
dipole hyperfme constant (A) for the ip 2 P\/2 state using different CC methods. 

have shown that unitary coupled cluster (UCC) method is capable of providing accurate estimates of 
wave lengths, transition probabilities, oscillator strengths, nuclear quadrupole moment, magnetic dipole 
and electric quadrupole hyperfme constants for relatively light atomic systems with a single valence 
electron. The calculated value of Q compared to others are closer to the experimental uncertainties 
than all the existing atomic and nuclear calculations, thereby demonstrating that RCC theory of atoms 
can yield accurate values of nuclear quadrupole moments. Such an inter-disciplinary approach involving 
atomic and nuclear physics adds a new dimension to this theory. 

Acknowledgments : One of the authors (CS) acknowledges the BRNS for project no. 
2002/37/12/BRNS. 
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